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Приближенный анализ многомерной 
конвективной математической модели 
кристаллизации металла 


Исследуется многомерная задача Стефана с учетом конвективного движения в жидкой фазе. Получена 
формула зависимости уравнения свободной границы от числа Рейнольдса. 


Постановка задачи 


Работа посвящена изучению процессов кристаллизации, когда распространение 
тепла связано не только с теплопроводностью, но и с конвективным переносом, реально 
присутствующим в жидкой фазе. 

Жидкая фаза рассматриваемого процесса заслуживает специального исследова- 
ния из-за априорной возможности существования поля скоростей, вызывающего интен- 
сивную теплопередачу путем конвекции, само же поле скоростей на самом деле воз- 
никает под действием объемных сил земного притяжения из-за разности плотностей 
жидкого шлака и металла, архимедовых сил из-за объемного расширения при значи- 
тельных температурных градиентах, а также в случае электрошлакового переплава, 
сильных электромагнитных сил, способных играть доминирующую роль и сообщать 
конвекции, наряду с силами тяжести, вынужденный, а не свободный характер. Усилен- 
ная циркуляция в расплавленной шлаковой ванночке еще в 1957 году была обнаруже- 
на в исследованиях академика Б.Е. Патона, академика Б.И. Медовара и их сотрудников, 
а дальнейшее экспериментальное изучение этих явлений, как отмечалось в докладе 
Ю.М. Каменского и др. на Международном симпозиуме по специальной электроме- 
таллургии (Киев, июнь 1972 г.), привело к выводу, что интенсивное перемешивание 
шлакового расплава имеет место по всему объему шлаковой ванночки и вызывается 
электромагнитными силами: под током средняя скорость перемешивания частиц шла- 
ка достигает 7 — 10 см/с, тогда как после отключения тока, несмотря на сохранение 
теплоотвода через стенки кристаллизатора, скорости падают до 1 см/с, то есть умень- 
шаются на порядок. Отмеченные величины скорости вынужденной конвекции срав- 
нимы с поперечными и продольными размерами шлаковой ванночки, следовательно, 
конвективная теплопередача в шлаковой ванночке значительно превосходит теплопе- 
редачу через теплопроводность. 

Основной темой данной работы является гидро- и термодинамические яв- 
ления в металлической, а не шлаковой ванночке, и поскольку здесь их закономернос- 
ти не менее сложны, а экспериментальные исследования, насколько известно, отсутст- 
вуют, мы исходили из предположения, что существует, по крайней мере, три зоны со 
своими характеристиками кинематических и температурных полей: 1) зона перегре- 
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того металла, поступающего в виде серии капель или сплошной струйки; 2) погранич- 
ный слой вдоль фронта кристаллизации с заданной на нем постоянной температурой; 
3) пограничный слой вдоль стенки кристаллизатора при наличии теплоотвода через 
нее. Последняя зона из-за теплоотвода через известную часть границы и наличия фа- 
зового перехода вдоль неизвестной границы придает модели температурного погранич- 
ного слоя большую теоретическую сложность, поэтому мы остановимся более подроб- 
но на результатах исследования первых двух зон. 

Для описания поля скоростей в зоне поступления перегретого металла исполь- 
зуется математическая модель затопленной струи вязкой жидкости, основанная на 
известном в теоретической гидродинамике точном решении нелинейной системы диф- 
ференциальных уравнений Навье-Стокса, а поскольку в затопленной струе перемеши- 
вание можно охарактеризовать так называемой свободной турбулентностью, харак- 
теризующейся одним числовым параметром — коэффициентом «кажущейся», или тур- 
булентной, вязкости У, , полученные формулы можно интегрировать как в ламинарном, 


так и в турбулентном приближениях. 
Пусть в пространстве А” задана область © ‚ граница которой состоит из двух 


замкнутых, связных, гладких, непересекающихся поверхностей Г’и Г` класса Н“*. 
Пусть далее Г, ‚ где 1 =[0,Т] — гладкие поверхности, лежащие внутри ©, и такие, что 


Г’ лежит внутри ограниченной области, границей которой является Г,, и Г, в мо- 
мент времени { разбивает © на две связные подобласти О! и О’. Двухфазная задача 
Стефана при наличии конвективных движений в жидкой фазе состоит в нахождении 
скорости жидкости У(х,г), давления р(х,/) , распределения температур и* (х,г) ии (х,0 
соответственно в жидкой и твердой фазах и свободной поверхности Г, по следую- 
щим условиям: 


+ а 
95 инф а2У2и+ (1) -0, (дер, 9% 650 122-51) -0, (дер, 
6 - т’ д Г 
’бЬЙ урн+Урсь д = УЗ. 0+ Ки*), УРо.д=0, (хде, 
0 Ве (1) 
+ + Е + > >. = 
и (>, 1) | =0 = А (х), и (х, 1) хеГАЛ- — В (х, 1), Г(х, 1) —0- |@/ 6.) й Г(х, 1) хЕГ+" Л, — 0, 
3 би ди* 
+ 
и (х,0 хеГ, =0, Ра а р соз(и,х;) +кс03(и,1) =0, хеГ,, 


1 1 
где 0; ={(х,): хе О е (0,Т)|,х= (мьх,,х,), У = (0/дх,0/дх,,0/дх,) ‚ ©. — области жид- 
кой и твердой фазы, на которые разбивает область © граница Г, ‚ причем &* =Г‹л*, 


п — нормаль к Г, ‚ направленная в сторону ©. Предполагается, что 
+ 3+ —- + .. 2+а 5+ + Е 
В—(х,) ЕН (Г х[0,7]) ›С(х) ЕН” (1), 4%) ЕН" (6%), 2) 
Хит) ССВ), +В) 2) >05 (д еГЁЖО,Т], 
где 0<В<а, ©х - области, на которые разбивает © граница раздела фаз Г, , а па- 
раметры а.,^,,к, Ве, =, — считаются положительными постоянными. Относительно 


функции /(и’) всегда можно считать, что ха) =0, так как общий случай сводится к 
и 


указанному после замены р(х,г) = р(х,6)+(1(1),х), 1(и)= (и) ГО. 
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Приближенное решение задачи (1) 


Для точек поверхности Г, введем координаты @ = (@,0,). Через х(@) Г, или 
через @ будем обозначать также соответствующие точки в А°. Известно, что Г, = 
= {х =х(@)+п(@). р(®,0} ‚где р(@,Г) — функция класса Н?*“”" (Г, х[о,Т]), р(@,0) =0, 


И + 
а п(@) — нормаль к Г, , направленная внутрь 62. [2]. 
Предложен метод решения задачи (1), состоящий в разложении решения в ряд 
по степеням малых чисел Рейнольдса Ке [1]: 


и" (х,;Ве) =иу(х)+ У (Ве (х,1), ро,вВе) = р) + Уве! Ре), 
ы я: (3) 
Г (х,5Ке) =И, (>) + У(во"и, (0), 1=1,2,3, р(@,6Ве) = У (Ве ри(.. 
Далее, если И ОИ выполненными ОЕ условия: 
У? А" (х)=0, хе 0$; 4" |. = В*(х,0), 4* В =0, С) =0, хеО!, 44 
|3 


УА- © =, 


У4* (|, 


Уд" (х)|. >2>0, 
0 
+ ня 
тогда в качестве нулевого приближения можно взять функции А“(х), т.е. из (х) = 
= 
= А*(х), хеО*. Причем из последнего условия (4) в предположении звездности по- 
верхностей Г“ следует, что поверхность Г, принадлежит классу С”, не имеющую 


самопересечений и расположенную относительно Г” и Г’ аналогично поверхности 
Г, в задаче (1) [2]. 


Построение первого приближения 


Далее, пусть р: =Оёх(О,Т), Г = Г*х[О, тр”. =Г,х[О,Т] . Тогда для первого при- 


ближения и’ (х,0), И (1) = (И, (х,6, И, (х,6), Г. (х,1)), р, (@,[) из условий (1) и разложе- 
ния (3) вытекает следующая задача: 


Урь(х) = УЗ (х, 0) + (из), (х,1) Е Бу, (5) 

УТ (х,1) = 0, (х,#) Е Б*; Г (х,1) нов, = Й(х,0 |. =0, (6) 

м (ур -аУи, (х,1) = Е" (к), (х,1) в Бу, иг (хо =0, (7) 

м (%1) |5 =0, [иго рко,д-+и 6). | =0, (8) 

и, (х,0 |; =0, [Учеб ре, +ия (6,0 | =0, (9) 

где я = - (И с-оу 509 при (х,г) еР+ иЕ (0) =0 При (х,0) ЕР». Задача (5), (6) фак- 


—_ 2+8,2+8 —_ „В 
тически изучена в [3] (теорема 3), причем РЕН о (5х); УРобдеН” 2 (5+). 


1+ ‚2+В 
Тогда Е .ЕН В, (5). 


46 «Искусственный интеллект» 272010 


Приближенный анализ многомерной конвективной математической модели... ЛМ 


Предполагаются также выполненными условия согласования до первого порядка 
включительно, которые вытекают как необходимые из предположения о существова- 
нии гладкого решения и формулируются аналогично [1, с. 363]. В частности доста- 
точно предположить: 


Уру(х)=0, Е (10) 


+ + + р 
так как аи) + Е, ща (х) |, (50). й = 
0 
2+а 1 
Далее, при заданном р (ЕН? 2) найдем функцию и (хе р)ЕН” 2 (2;} 


как единственное решение задачи (7), (8) [1]. Построим затем оператор М, следую- 


щим образом: 


М. р, Е 


ди 


Вы 
[& > (х(@), 5 р,)- К. В 


( х(6), 1; р, › х(@)>ЕГ., который действует из 
т 


+ (@+2) 
и 
2: 


2+а,^"% ое и 
Н ’? ИЯ в т а (2. Для решения и*(х,1; р’) справедливы оценки 


(&+2) 


(а) 
<е(Е* +| 
( 1 ря и 


), где с — некоторая постоянная [1, с. 364]. 


2+а 


0 раке 
Рассмотрим функцию р (о, е А (Р,) и построим соответствующее реше- 


й О 
ние и! =и’(х,6; р,). Тогда получим 


о |(@+2) (&+2) 


+ 


+ 
и си 


0 


<ер - р, 


Гт 


Отсюда, с учетом того, что Мр-м Е ы -| Е д(и; —и;) й д(и —и, ) + › сле- 
1771 Я 


ие 108 0 


0 
‚ где с=с(к, +К_)Т/К. 


(+2) (@+2) 


| 


М.р,-М, р, 


0 
дует 


Гт Г 


Итак, оператор М, сжимающий, если выполняется условие 
СЕТА 1, (11) 
Это неравенство всегда выполнимо, например, при малых значениях Т. Поэто- 
+а 
И Г М Бе 
му оператор М, имеет неподвижную точку в Н ( ‚ ‚ те. Ир = ру. 
Лемма 1. Пусть выполнимо условие (11). Тогда оператор М,, действующий из 
2+а 2+а 
2+а,—— 2+а—— 
то [в в” Е. ‚ имеет там неподвижную точку. 


Теорема 1. Пусть выполнены условия (2), (4) и (10), тогда существует единст- 


2+а 


венное решение и"(х,г), И (х,0, р; (©,0) задачи (5) — (9), причем и*(х,1) ый Чо (р, ), 


2+в,2+В. 


7(х,еЕН (2) ир (®0Е Е (Г) , 
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Далее рассмотрим второе приближение 7 (х,1) ; из (х, о, Р> (©, г) задачи (1) для 


малых чисел Ве. Имеем: 


УР! (1) = У2У, + Г м *, 
+ 


ди> 
7.0509 =0; р 


[6] 


ВУ РНЕ 
аз У^и> 


— р 
=^ (х,#) в 


р 


+ = 
Т› 4 с =0, (12) 


[Умер (дни 0,0) |= ЛЕбь)), хеГо, ивбы д =0, 


+ 
ди> 


[4 —^ 


би> 


+ ди 


[611 


др 
е 2 
= 5 + 1,1), хеГ 


0? 


где Ес) =-(ИУ мг - (У) при (хбеЕР: и Е (х,0 =0 при (х/еЕВ?, ое 


Е 5 Есцо Ор.) 
И 


Рассмотрим оператор М. : 


1, ды 
М.р, = Н 2 (х(@),;р, 
о 


ди 


1’ (х(@) +тп(®)р,(@,1)) 


Ат” 


‚ хЕГ.. 


нк Я (щоувр.)- побои х(@о)Е Гу. 


Лемма 2. Пусть выполнены условия леммы 1. Тогда оператор М, имеет в 


24а, 2% 


Н - (2) неподвижную точку. 


Для задачи (12) условие согласования первого порядка выглядит следующим 


образом: 


УР (х,0) =0,х(@) Е Го. 


(13) 


Теорема 2. Пусть выполнены условия теоремы 1 и пусть выполнено условие (13). 


ЕО 
Тогда существует единственное решение задачи (12), причем и:(х/)еН `? (2; 


5 +в,2+В. й 2+ 
И<.)еЕН” 2 (5, } р. (@,) ЕН 


"2 (Р,), Ур о.деН” 


2(В:). 


О влиянии конвекции на фронт кристаллизации 


Основной целью дальнейшего исследования задачи (1) является изучение гидро- 
динамических явлений в жидкой фазе. Справедливо утверждение. 

Теорема 3. Пусть выполнены условия теоремы 2. Тогда при малых числах Рей- 
нольдса Ве и достаточно малых значениях { справедлива формула: 


Г;:х=х(@)-Веп 


из (*(6),0) 


МА (х(6)) 


2:1), - Де), 
У4(х(6)) 


(14) 


+0 (ве? |) Е Г, 


где и (х, Виз (х,г) — решения соответственно задач (5) — (9) и (12). 
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Последняя формула позволяет исследовать свободную поверхность Г, в зависи- 


мости от чисел Рейнольдса и проследить, как конвекция влияет на процесс кристал- 
лизации, что проверить в лабораторных условиях практически невозможно. 
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О.С. Мшенко 

Приблизний аналВ багатовимрно! конвективно! математично! модел! кристал зацй металу 
Дослджуеться багатом!рна задача Стефана з урахуванням конвективних рухв у рдиннй фаз1. Доведена 
формула залежност! р1вняння вльно! границ в1д числа Рейнольдса. 


А.5. Мтепко 
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